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Abstract. In this paper we deal with the noteworthy Sylvester’s determinantal identity and some 
of its generalizations. We report the formulae due to Yakovlev, to Gasca, Lopez-Carmona, Ramirez, 
to Beckermann, Gasca, Muhlbach, and to Mulders in a unified formulation which allows to understand 
them better and to compare them. Then, we propose a different generalization of Sylvester’s classical 
formula. This new generalization expresses the determinant of a matrix in relation with the determinant 
of the bordered matrices obtained adding more than one row and one column to the original matrix. 
Sylvester’s identity is recovered as a particular case. 

1 Introduction 

Sylvester’s determinantal identity |22| is a well-known identity in matrix analysis which 
expresses a determinant composed of bordering determinants in terms of the original 
one. It has been extensively studied, both in the algebraic and in the combinatorial 
context, and it is usually written in the form used by Bareiss [3|. 

In the last years some generalizions of this identity have been proposed, with useful 
applications mainly in interpolation and extrapolation problems ma Eg. One of the 
first generalizations seems to have been proposed by Yakovlev mm who used it in 
the control theory domain. In m, Gasca, Lopez-Carmona, and Ramirez proved a 
generalization of Sylvester’s classical determinantal identity; an alternative proof based 
on elimination techniques was suggested in Eo], while a short and simpler one than 
the earlier attempts was given in [6] . More recently, Mulders m formulated another 
generalization related to fraction free Gaussian elimination algorithms. Moreover, in 
order to simplify the reading when dealing with minors identities, Leclerc m proposed 
a unified notation to encode many algebraic relations between minors of a matrix, 
including Sylvester’s determinantal formula and the Mtihlbach-Gasca generalization 
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These generalizations (and sometimes only a rewriting of Sylvester’s formula as 
Yakovlev’s) are scattered in the literature. Moreover, they used different notations 
which makes difficult to understand them as a whole and to see the connections be¬ 
tween them. Therefore, one of the important aims of this paper is to present these 
identities in a unified notation. Also, for helping the reader to better understand these 
generalizations, they are illustrated by several examples. 

Then, for showing the potentiality of our unified approach, a new generalization of 
Sylvester’s identity will be given. Sylvester’s determinantal identity can be recovered as 
a particular case of it. Such identities may have a potential impact for the development 
of new recursive algorithms for interpolation and extrapolation. 

In Section 2, the classical Sylvester’s determinantal identity is reminded, and some 
particular cases are discussed. Section 3 is devoted to the generalizations, mentioned 
above, of the classical identity. These identities are presented in a unified formulation 
which allows to understand them better and to compare them. In Section 4, the new 
generalization is proposed. 


2 The classical Sylvester’s determinantal identity 

To introduce the classical Sylvester’s determinantal identity we consider a square matrix 
M = ( a,ij ), of order n over IK (a commutative field). We denote its determinant by det M 
and, for an integer t, 0 < t < n — 1, and a given couple of integers i,j with t < i,j < n, 
we define the determinant as 


On • 

• an 

°i j 

021 • 

■ a2t 

02 j 

an • 

att 

a tj 

Mil 

^it 

Clij 


for 0 < t < n 


1 and a[°- = a^, 


(1) 


that is the determinant of order (t + 1) obtained from the matrix M by extending its 
leading principal submatrix of order t with the row i and the column j of M. Given 
these assumptions we can easily write down, in the following Theorem, the Sylvester’s 
identity. 


Theorem 1 (Sylvester’s identity) Let M be a square matrix of order n and t an 
integer, 0 < t < n — 1. Then, the following identity holds 


det M ■ 


At- 1 ) 

h,t 


1 n—t—1 


a 


(t) 

't+i.t+i 


a 


a. 


A) 

1 


w 

t+l,n 


,(*) 

l n,n 


( 2 ) 


with ag q 1J = 1. 


For different proofs of this Theorem, we refer to mmm- We note also m , where 
a combinatorial proof, together with several non-commutative extensions, has been 
presented. 
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Some particular cases of this Sylvester’s identity are of interest: 


Case 1: For t = 1 and an / 0 we obtain the Chid pivotal condensation method m 
which states 


det M 


det B 
(on) n_2 


where B = (&*•/), is the square matrix of order (n — 1) with entries 


f>ij 


= a 


(i) 

i+lj+l 


ail °i,j+i 
Oi+1,1 Oi+l,j+l 


for = 1 ,... ,n — 1 . 

Case 2: For t = n — 2, we consider M partitioned like 


/ a h T e \ 

M= c D f (3) 

\« h T l/ 

where a , e, g, l are scalars, b, c, f , h are (n — 2)-length column vectors, and D is 
a square matrix of order (n — 2). It is easy to apply the Sylvester’s identity to 
such a matrix. In fact, by interchanging rows and columns in order to put D in 
the upper left corner, leaves the sign of determinant unaltered, so 


det M = 


and applying ([2]) we have 


det M • det D = 


D c 

b T a 

D c 

h T g 


By setting 
A' = 


a b 1 

c D 


B' = 


b T 

D 


D c f 

b T a e 
h T g l 


D f 

b T e 

D f 

h T l 


, C' = 


a b T 

c D 

c D 

g h r 


c D 

g h T 


, D' = 


b T e 

D f 

D f 

h 7 l 


D f 

h T l 


we can write the very simple rule 

det M det D = det A' det D' — det B' det C', 


( 4 ) 


( 5 ) 


frequently used to obtain recursive algorithms in sequence transformations (see, 

for instance, Ezusum). 
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3 Generalizations of Sylvester’s identity 

Several authors have deepened the main property of classical Sylvester’s identity, which 
is very useful in several domains. Some of these authors have generalized the classical 
identity obtaining interesting new formulae. Yakovlev [23; 2M proved a generalization 
useful in control theory. After, Beckermann, Gasca, Miihlbach et al. 13 EH Eo] devoted 
some papers to generalizations and applications of Sylvester’s identity, while Mulders 
[21 j proposed a particular generalization of practical use in linear programming and 
related topics. In this section, we present a brief overview of these results. 

In order to have a common working framework, we first introduce here some general 
notations and definitions related to a n x m matrix M over IK. 

An index list I of length k < n is a fc-tuple (with possible repetition) of integers 
from {1, 2,..., n}, that is / = (i\,... , i &), with 1 < ii < n for l = 1,... k. If we impose 
that all the elements are different, we have an index list without repetition. So, in this 
case, card(I) = k. If, in addition, a < /3 implies i a < ip, then / is called an ordered 
index list. The usual symbols €,C,C,U,D,\, used for the sets will be used for the 
ordered index lists and the result will be again an ordered index list. For any n £ N + 
(the set of all positive integer), we define the ordered index list N n = (1,2,..., n). 

In the sequel, if nothing is stated to the contrary, all index lists are assumed to be 
ordered and without repetition. 

Let now I = (i\,... , i a ) C N n and J = (j\,... ,jp) C N m two ordered index lists. 
We denote the ax/3 submatrix, extracted from M, with rows labeled by i £ I and 
columns labeled by j £ J as 


M (T) 

= M l 1 } ’' 

• j *aA_ 

^ a *u'l 

a hjp ^ 

W 

yi! • 

■>W - 


a iajp / 




\ a i a jl 


( 6 ) 


Let us mention that, in some papers, the indices for rows and columns are inverted. 
However, the notations we are using, seem to be the most common ones; see, for 
example, [1J. 


If a = j3, n = m, we denote the corresponding determinant as 








M 

y 

= M 

h, 

• • 5 ^ a 

— 


j 


ji, 

• • ? ja 


notations 

(t) 

a) ~ = 

M 

!;••• 

t, i 






tj 


-Hu 1 


i 




'I'OiJoc 


, for t < i,j < n and 0 < t < n — 1 . 


3.1 Generalization of Yakovlev 

In [22, 21], Yakovlev presented some determinant identities used to simplify several 
procedures in control theory. In particular, in [21] he was the first to give a simple 
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generalization of Sylvester’s identity, which can be viwed as a renaming of the row and 
column indices. 

Assume that M = (a,;j) is a square matrix of order n. Let 0 < t < n— 1, and consider 
the two ordered index lists I = C N n , J = (j\.... .jt) C N n , and the 

complementary ordered index lists I' = {i' t+1 ,... ,i' n ) C N n , J' = (j' t+ ,,... ,j' n ) C N n 
(that is / U I' = J U J' = N n ). 

Starting from the well-known Laplace expansion formula for the determinant of a 
squared matrix, as proposed in EH, the author proved a determinantal identity and 
use it in the solution of certain applied problems as regulator synthesis for system of 
hight order. This identity can be transformed in the following formula 


det M ■ 



h, ■ ■ ■ ,it 
Jl, * * ’ 5 jt 


n—t—1 


n 

e (- 1 )" n m 

P{ptt + l, — ,Ot n ) f)=t+l 


*1 5 ' ' ' , ^OLp 

Jl, * * * , jt 1 j ft 


( 7 ) 


where P(at+i, ■ ■ ■, a n ) represent the set of all permutations of (t + 1,... , n) and [i is 
the number of inversions needed to pass from (t + 1,..., n) to a certain permutation 
{ptt+ 1 1 ■ • • i ®n)• 

Let us now take I = J = (1,..., t). Thus, I' = J' = (t + 1,..., n). Formula (|T|) 
becomes 


det AT - (m 

'1 

1. t 

\ n—t—1 

) = s 

n 

(-i r n M 

1, • 
1, • 

. , OL[3 

..t.3 

\ 

_ ? J _ 

P(at+ i,...,c 

P=t+1 

5 

) ) r 


This expression is one of the different ways to formulate the Sylvester’s determi¬ 
nantal identity. In fact 

1, . . . , t, Otp 

t,f3 

and, by Leibniz formula, 


M 


(t) 

= a 4,f>’ 


e (- 1 )" n 1 

P((*t+l,---,Oln) P—t+1 


l {t) 


is exactly the determinant of order n — t on the right hand side of ((!ZJ). 

In the following examples we consider P(ctt+i,... ,a n ) as an ordered set of (n — t)! 
elements, and we denote by /a = (/xi,... ,Ht n _t)\ ) the vector whose elements are the 
corresponding number of inversions. 

Example 1: 

Let n = 6 the order of M, and set t = 4. We consider the ordered index lists I = 
= (1,3,5,6), J = (ji,...,j 4 ) = (1,2,4,6); then I' = (i' 5 ,i' 6 ) = (2,4), 
J' = C?5, Je) = ( 3 , 5 )- We have 


det M ■ M 


E (- l ) M II M 

P(a 5 ,o 6 ) 


3=5 


i a 

1 •> L c 


j, 
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Since P(a^, oq) = {(5,6), (6,5)}, and p = (0,1), we obtain 
det M ■ M 


= M 


I, 2 

J, 3 


■M 


IA 

J, 5 


- M 


IA 

J, 3 


• M 


IA 

J, 5 


With the same values, but with I = J = (1,..., 4), I' = J' = (5, 6), P{a 7) ,aA = 
{(5,6), (6, 5)}, and /x = (0,1) we obtain 


det M ■ M 


= M 


I, 5 

J, 5 


M 


1,6 

J, 6 


- M 


1,6 

J, 5 


M 


IA 

J, 6 


that is 


det M ■ M 


X- 

■ A 


1,- 

■A_ 



which is the classical Sylvester’s identity. 


( 4 ) 

<5 

( 4 ) 

a 6,5 


( 4 ) 

°5,6 

( 4 ) 

6,6 


a, 


Example 2: 

Let now n = 8 and t = 5, and consider the ordered index lists I = = 

(2,3,5,6,8), J = (Ji,... ,jA = (2,4, 5,6,7), and the complementary /' = = 

(1,4,7) J' = (j' 6 ,f 7 ,f 8 ) = (1,3,8). We have 



Since P(a 6 ,a 7 ,a 8 ) = {(8, 7,6), (8,6, 7), (7, 8, 6), (7, 6, 8), (6, 7, 8), (6,8, 7)}, and = 
(3, 2,2,1,0,1), we have 


det M • (^M 


) =_„ 

IA 
JA 

■ M 

IA 

J, 3 

• M 

IA 
J, 8 

+ M 

1,7 

JA 

■ M 

IA 

J, 3 

• M 

,00 



M 

IA 
JA 

■M 

I, 7 

J, 3 

■ M 

IA 
J, 8 

- M 

IA 
JA 

• M 

IA 

J, 3 

• M 

00 1 
S'. 



M 

IA 
JA 

■M 

CO 1 

S', 

■ M 

I, 7 

J, 8 

- M 

IA 
JA 

• M 

I, 7 

J, 3 

■ M 

’i-H 

,00 


3.2 Generalization of Gasca, Lopez-Carmona, and Ramirez 

In ! 15i . Gasca, Lopez-Carmona, and Ramirez proved a first generalization of Sylvester’s 
determinantal identity; this generalized identity has been applied to the derivation of 
a recurrence interpolation formula for the solution of a general interpolation problem. 
Assume that M is a square matrix of order n, and that n = t + q, for some positive 
integers t, q. Given a set of q ordered index lists Jk = (j\,... ,jt+i) C A T n , k = 1,..., q, 
with 

J card( Jfc) =7 + 1, k = 1,..., q, 

\ card (n J k + 1 ) = t, k = l,...,q-l, 

= u j k 

k =1 

Sjc t/fe n Jk -\-1 (#1,..., ), k 1? • • • j Q l* 


we set 
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Let B the matrix with elements 
b ik = M 

and let be the element of Jk such that 


1 , 2 ,t,t + i 
Jk 


l < i,k < q 


•ihu ^ Jk —1 Jk Jk 1 ■ 1 2 , . . . , q 


jft , e Ji — Ji — J\ — <51, 


k = l. 


With these notations the authors give the following generalization of Sylvester’s identity 


9-1 


det B = c • det M • M 


fc=i 


1 , - - -, t 

S k 


( 8 ) 


where c is a sign factor which does not depend on the element a t] of M, but only 
on the set of Jk according to 


c = 


with ^ = q(q - 1)/2 + J2(j k hk - h k ) 


k= 1 


if card(J^) < t + q 
if card(J^) = t + q 


<? 

with J^ = U Jk- 

k=1 


When Jfc = (1,..., t, t + k), for k = 1,2,... ,q, we have, for all k, Sk 
and the identity (| 8 |) gives 


( 1 ,■■■,*), 


det B 


det M ■ 



q -1 


which is exactly ([ 2 ]), when n = t + q. 

An alternative proof of this identity was suggested by Gasca and Miihlbach in [2U]; 
this proof is based on the Laplacian expansion and on an elimination strategy. Roughly 
speaking, when a determinant is simplified by performing elementary operations leaving 
its value unchanged, then elimination strategies are applied, for example Gauss or 
Neville elimination. 


3.2.1 Examples 
Example 1: 

Let M be a nonsingular square matrix of order n = 5, choose t = 2, q = 3, and 
consider the ordered index lists J\ = (1,3,4), = (1,4, 5), J 3 = (2,4,5) and j( 3) = 
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(1, 2,3,4, 5) = IV5. We have Si = (1,4) and S 2 = (4, 5). Let 



( M 

1.2.3 

1.3.4 

M 

1,2,3 

1,4,5 

M 

1,2,3 

2,4,5 

\ 

B = 

M 

1.2.4 

1.3.4 

M 

"1,2,4 

1,4,5 

M 

1.2.4 

2.4.5 



\ M 

'1,2,5 

1,3,4 

M 

"1,2,5 

1,4,5 

M 

'1,2,5 

2,4,5 

/ 


Since card(J^) = n, for computing the sign factor we hrst have to determine jfc for 
k = 1, ■ ■ ■, 3. For k = 1, j' hl = 3, hi = 2, for k = 2, = 5, h 2 = 3, for k = 3, j| 3 = 2, 

/13 = 1. So /r = 3 + 4 = 7 and c = —1. 


Finally we have 


det-B 


det M ■ M 


1,2 

1,4 


•M 


1,2 

4,5 


Example 2: Let again n = 5, £ = 2, q = 3, but consider the ordered index lists 
Ji = (1,2,3), J 2 = (2,3,4), J 3 = (1,2,4) with card(j( 3 ) = (1,2, 3,4)) <5. So c = 0, 
and the determinant of the matrix 



( M 

1,2,3- 

1,2,3 

M 

1,2,3" 

2,3,4 

M 

1.2.3 

1.2.4 

\ 

B = 

M 

1,2,4" 

1,2,3 

M 

1.2.4 

2.3.4 

M 

1,2,4" 

1,2,4 



V M 

1,2,5" 

1,2,3 

M 

1,2,5" 

2,3,4 

M 

1,2,5" 

1,2,4 

/ 


is equal to zero. 


3.3 Generalization of Beckermann, Gasca and Miihlbach 

Sylvester’s classical identity can be interpreted as an extension of Leibniz’s definition 
of a matrix determinant, known as the Muir’s law of extensible minors CD- Muir hrst 
stated the Law of Complementaries, which was been already known to Cayley in 1878 
and reads as 

To every general theorem which takes the form of an identical relation between a 
number of the minors of a determinant or between the determinant itself and a number 
of its minors, there corresponds another theorem derivable from the former by merely 
substituting for every minor its cofactor in the determinant, and then multiplying any 
term by such a power of the determinant that will make the terms of the same degree. 

Brualdi and Schneider [12] gave a formal treatment of determinantal identities of 
the minors of a matrix and then provided a careful exposition of the Law of Exten¬ 
sible Minors and of the Law of Complementaries, two methods for obtaining from a 
given determinantal identity another determinantal identity, as mentioned earlier. In 
m , Miihlbach using the relation between Muir’s Law of Extensible Minors, Sylvester’s 
identity and the Schur complement, presented a new principle for extending determi¬ 
nantal identities which generalizes Muir’s Law. Miihlbach’s proof makes use of general 
elimination strategies and of generalized Schur complements. As applications of this 
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technique, he derived a generalization of Sylvester’s identity, which is the same as the 
one proposed by Gasca and Miihlbach in [20|. Later, Beckermann and Mtihlbach [6] 
gave an approach, shorter and conceptually simpler than the earlier attempts mm, 
for a general determinantal identity of Sylvester’s type. In the sequel we report these 
last results. 

We consider a n x m matrix M with rows and columns numbered in the usual way. 
Let I = C N n and J = (ji, ■ ■ •, jp) C N m two ordered index lists, and 

consider also the row vectors zi,. .. ,z q , where zy, = 

We border the matrix M , extracted from M, by the q row vectors z' k = 
( z k,j\)\=i,...,p , k = 1,q, extracted from the z*,, that is we consider 

/ a h,ji ''' a h,jp ^ 


M 


*1) ■ ■ • da | 1) • • • j Q 


31, 


,3/3 




l la,3/3 


Z fj 1 


Z l,30 


V z qJi ''' z q,jp / 

Let q £ N + be fixed, and consider the index lists 

II, • • • , Iq C Nn and J \. . . . . Jq d N m 

such that 

card(Jfc) = card(Jfe) + 1 for k = 1, ..., q. 
We set, for k = 1,..., q, 

k k 

/( fc ) :=f)li and J (fc) := (J .J 3 


2 — 1 


j = l 


The authors give the following theorem. 


Theorem 2 Let 


Iq Q I^ such that card(/o) = card(J < - 9 ^) — q. 

Then, for any matrix N n x N m matrix M, there exists an element c £ IK depending only 
on M suc -h that, from all row vectors zi,... ,z q £ and so for the vectors 

z'k = ( z k,j\)j x £j(i)> k = 1,... ,q, extracted from them, 


det B = c • M 

where B is the matrix with elements 

Ii\j 
Ji 


I 0 \l,...,q 


( 9 ) 


bij = M 


1 < i,j < q 
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that is the matrix M 



bordered with the row vector z'■ 




Two corollaries in [5] are also of interest. 

Corollary 1 If card(/^ fc ^) > car d(J^) - k, for some k € {2,... , q}, then (G|) holds 
with c = 0. 

Corollary 2 If card(J^) < q, then, for all choices of vectors z\,...,z q , we have 
det B = 0. 


In some particular cases, the constant c of Theorem [2] can be computed explicitly. 
This is the case of Sylvester’s classical determinantal identity. In fact, we consider 
a nonsingular matrix M of order n, and we choose two integer t,q £ N + , such that 
t + q = n. Let 


Jk — 

(i, ■ ■ 

-,t,t + k), 

for k = 1,.. 

• id 

h = 

(1, ■ ■ 

-,t), 

for k = 1,.. 

■ ,q 


/ 0 = /(9) = (l,...,t), J(9) =N n = (!,..., n), 


z' k = M 


t + k 
Jk 


for k = 1, • • •, q- 


We consider the matrix B with elements 


hj = M 


!,•••, t\i 

l,...,t,t+j 


From (l9l) we have 


det B = c • M 


l,...,t + q 
l,...,t + q 


1 <i,j <q- 


= c ■ det M 


( 10 ) 


where c depends only on M 


h 

N m 


and can be evaluated by choosing the particular 


row vectors = (5x,t+k)xeN n for k = 1,..., q, where 5 is the usual Kronecker symbol. 
With this choice we have 

det B = (m I" 1 ’ - "’* 


1 


while the right hand side becomes 


* * * ? 1 


So, assuming that M 


1 

1 


+ 0, 


c= M 


1 

t 


q -1 
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and from (1101) it follows 


det B = det M ■ yM ^ 

which is the classical Sylvester’s identity ([ 2 ]). 

In a similar way [ 6 ], is also possible to obtain the identities of Schweins and Monge 

HI- 

The relation Q can be very useful for obtaining other determinantal identities, 
since c depends on M but not on the rows corresponding to indexes from N n \ Iq. 
As an application of this general determinantal identity, in [ 6 ] old and new recurrence 
relations for the E-transforms [7J are derived. In the same paper, a more general 
Theorem, based on a chain of index lists (that is lists Ik C N n and Jk C N m such that, 
for all k = 1 ,..., q , card( J^) = card(/( fc )) + k) is also provided. 



3.3.1 Examples 
Example 1: 

Let M a n x m matrix, with n = 6 and m = 7. Set q = 3, and consider the ordered 
index lists/i = (2,3,4 ),/ 2 = (2,4 ),/ 3 = (1,2,4) and = (2,3,4, 5), J 2 = (2,3,4), J 3 = 
(2,3,4, 7). We have I® = (2,4) and J ® = (2,3,4, 5,7). Choose any zi,Z 2 ,z 3 , and 
/ 0 = /( 3 ). Let 


/ 


B = 


\ 


M 

M 

M 


'2,3,4| 1 
2, 3,4, 5 
'2,4 11 

2,3,4 
1,2,4|1 
2, 3,4, 7 


M 

M 

M 


2 . 3.4 | 2 ' 

2,3,4, 5 

2.4 | 2 

2.3.4 

1.2.4 | 2 

2,3,4, 7 


M 

M 

M 


2.3.4 | 3' 

2,3,4, 5 

2. 4 | 3' 

2.3.4 

1.2.4 | 3 

2,3,4, 7 


\ 


'2,4 11,2,3' 
[2,3,4, 5, 7 
priori) such that © holds. 


where det B / 0. If M 


7 ^ 0 , there exists c (which cannot be computed a 


Example 2: Let again n = 6 and m = 7, but consider q = 5 and the ordered index 
lists Ii = (1,4), I 2 = (3,4), I 3 = (2,4 ),/ 4 = (4, 5 ),/ 5 = (4,6) and Ji = (1,2,4), J 2 = 
(1,4,7), J 3 = (1,2,7), = (1,4,7), J 5 = (1,2,7). We have I& = (4) and j( 5 ) = 

(1,2,4, 7). Let 



( M 

'1,4 11 
1,2,4 

M 

'1,4 1 2 ' 
1,2,4 

M 

"l, 4 1 3 
1,2,4 

M 

"l,4 4" 

1,2,4 

M 

1.4 | 5' 

1.2.4 

\ 

B = 

M 

3,4 | 1 
1,4,7 

M 

3,4 | 2 
1,4,7 

M 

3,4 | 3 
1,4,7 

M 

3,4 | 4 
1,4,7 

M 

3,4 | 5 
1,4,7 



\ M 

2,4 | 1 
1,2,7 

M 

2,4 1 2 
1,2,7 

M 

2,41 3 
1,2,7 

M 

2,4 | 4 
1,2,7 

M 

2,41 5 
1,2,7 

/ 


Due to CorollaryEl since card(jl 5 ') < 5, then, for any z 3 ,..., Z 5 , det B is equal to zero. 
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3.4 Generalization of Mulders 


Another generalization of Sylvester’s identity is due to Mulders [ 23 ]. The author follows 
the idea of Bareiss mm, where the Sylvester’s identity is used to prove that certain 
Gaussian elimination algorithms, transforming a matrix into an upper-triangular form, 
are fraction-free. Fraction-free algorithms also have applications in the computation of 
matrix rational approxinrants, matrix GCDs, and generalized Richardson extrapolation 
processes [5j. They are used for controlling, in exact arithmetic, the growth of inter¬ 
mediate results. The generalization of Sylvester’s identity, due to Mulders, was used to 
prove that also certain random Gaussian elimination algorithms are fraction-free. The 
fraction-free random Gaussian elimination algorithm, obtained in that way, has been 
used in the simplex method (to solve linear programming problems) and for finding a 
solution of Ax = b, x > 0, with A a m x n matrix. 

Given a n x m matrix M, we consider the index list / = (i\,... ,it) C N n (without 
repetition but not necessarily ordered) and another index list J = (j i,..., jt) C N m , 
(with possible repetition and not necessarily ordered). The two lists have the same 
length t , 0 < t < min(m, n). Notice that it is possible to have t = 0, so an empty list. 
With these assumptions, we can define an equivalence class [(zi, Ji),..., ( it,jt)\ of all 
possible permutations of the t pairs and define the determinant 




hi ■ ■ ■ At 
ji, • • •, Jt 


with 



1 when t 


determinant is equal to zero. 


0. Note that, when J has repeated elements, the 


We also define the operation on the pairs defining the above determinant as 


[(«i,vi), ■ • •, (u r ,v r )]. 

I and U = (ui,..., u r ) are index lists without repetition but not necessarily ordered, J 
and V = (ui,..., v r ) are index lists not necessarily ordered but with possible repetitions. 
This operation produces a new class of pairs by adding to the first class all the pairs 
( Uk,Vk ) of the second class when Uk 0 /, and by replacing the pair (ig,je) of the first 
class by the pair (uk,Vk), when Uk = ig for some ig € /. This operation is performed 
for k = 1,..., r, where r is an integer. This operation is associative and commutative if 
/ n U = 0. For example, the result of [(2,1), (3,3), (1, 5), (5, 3)] [(1,1), (2, 3), (4,4)], 

is [(2, 3), (3, 3), (1,1), (5, 3), (4,4)], since the pair (4,4) has been added, and the pairs 
(1,1), (2, 3) have replaced the pairs (1,5), (2,1), respectively. 

We consider now two integers i and j, l<i<n,l<j<m and the determinant 

a \{ ' (il) 


Due to the operation previously defined, the result is the following: 


• When i 0 /, then we simply extend the matrix with the zth row and jth column 
of M, and we compute the determinant of the new square matrix of order (t+ 1), 
that is 




H i • • • i H i ^ 

jl 5 • • • 1 jtl j 
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• When it exists 1 < k < t such that = i, we replace, in cm* 1 ’- 7 i)>—the 
pair ( ik,jk ) by the pair (i,j), that is we take the following determinant of a new 
matrix of order t 




i \; • • • i ik — i }••• jit 
Ji ? • • • ? .1 ? ■ ■ •: jt 


Let us consider now the particular case where = k,jj~ = k, for all k = 1,...,!, 
that is / = J = AL When t < i,j < min(n, m), we recover the usual definition 
afj given in ©, but not in the other cases. For simplicity, in the sequel, we set 
af] = a^ 1,1 )’—w e have 


If t < i,j 
If j <t <i 
< If i,j<t,i^j 
If i = j <t 
, If i < t < j 


of order t + 1 
of order t + 1 
dp- of order t 
dpj of order t 
dp] of order t 

L iJ 



5=0 
a$ = 0 

dp = O-K 1 ’!)’—’(*’*)] 

df\ is a new determinant. 


By using this particular choice for I and J, Mulders formulated the following gen¬ 
eralized Sylvester’s determinantal identity that allows to consider the determinantal 
elements afl when i and/or j are < t. In the sequel, we state a more generalized 
identity. 


Theorem 3 For 0 < t < min(n, m), 0 < p, q < t and 1 < s < min(n — p,m — q), the 
following identity holds 


a [(l,l),...,(t,t)]-f-[(p+l,q+l),...,(p+s,g , +s)] 


h,t 


s —1 


~(t) 

a p+l,q+l 

~{t) 

a p+l,q+s 

~(t) 

a p+s,q+1 

• Cl v ' 

p-\-s,q-\-s 


( 12 ) 


For m = n,p = q = t and s = n — t we obtain the classical Sylvester’s identity 


3.4.1 Examples 
Example 1: 

Let M a n x m matrix, with n = 7 and m = 8. Set t = 5, and I = J = N§. Assume p = 
3, q = 4, and s = 3, so in this case p+s > t. Then we have a [( 1 ’ 1 )>--->(5,5)]-(-[(4,5),(5,6),(6,7)] _ 

M (d t’ ^ and a!p\ = M ''' ’ ^ . The right hand side of (fT2l) is the following 


'1,2, 3,4, 5,6' 

and = M 

1, • • ■,5 

1,2, 3,5, 6, 7 


1,...,5 


determinant of order 3 


“(5) ~(5) ~(5) 
x 4,5 


*4,6 


*47 


a 


(5) ~(5) ~(5) 


5 5 ^5 6 ^5 'j 

(5) ~(5) ~(5) 

6,5 u 6,6 u 6,7 


ak 


a, 
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that is, thanks to the operation defined, 


M 


0 

M 

'1,2,3,4,5 

M 

"1,2,3,4,5" 

1,2,3,6,5 

1,2, 3, 7, 5 

1,... ,5 

M 

'1,2,3,4,5 

M 

"1,2,3,4,5" 

_1,...,5 

1,2, 3,4, 6 

1,2, 3,4, 7 


0 


M 


1 , - • - ,6 


M 


1,2,3,4, 5, 6 
1,2,3,4, 5, 7 


and (fl2|) holds. 

Example 2: 

Let again M a 7 x 8 matrix. Set now t = 6, and I = J = N%. Assume p = 2, q = 3, and 
s = 3. In this case p + s < t, so we have a [(hi),-,(6,6)]e-[(3,4),(4,5),(5,6)] = M 


1.2.3.4, 5, 6 

1.2.4, 5,6, 6 


I,---, 6 

1, - - - ,6j 

presence of two equal columns in the first determinant. The right hand side is now 


and 41 = M 


It follows that the left hand side is equal to zero for the 


S (6) 

u 3,4 

s (6) 

a 3,5 

~(6) 

~(6) 

a 4,4 

a 4,5 

~(6) 

~(6) 

a 5,4 

a 5,5 


a (6) 

u 3,6 

h (6) 

a 4,6 

h (6) 

a 5,6 



M 

0 

'I,---,6' 

— 

1,... ,6 



0 


M 


0 

0 

1 ,... ,6 
1 , • • - ,6 


= 0 , 


and (fl2l) holds. 


4 A new generalization 


We present now a different generalization of Sylvester’s identity. Let M be a matrix 
of dimension n x m. Chosen two indices t and s, 0 < t < r, with r = min (n,m), and 
1 < s < min(n — t,m — t), in such a way that it exists q £ N which satisfies r = t + q s, 
we can build a sequence of submatrices of M, 


M k = M 


(I, ■ ■ ■ ,t,t + 1,... ,t + ks\ 
\1,... ,t,t 1,... ,t ksj 


for 0 < k < q 


(13) 


(t—1) (r— 1) 

where, by definition, detMo = a) t ' and det M q = a^ r ■ We avoid to indicate 
the dependence from t and s since they are fixed. So, for the same reason, we set 
m k = t + ks, for k = 0,... ,q — 1, and we consider the matrix B k of order s, with 
elements 



1,.. .,m k ,m k + i 
1,.. .,m k ,m k +j 


_ (rn k ) 

u m k +i,m k +j > 


1 < i,j < s. 


(14) 


P 

We remind that, by definition, when a > /3, the product JJ © is empty, and its 

i=a 

value is set to 1, and that, when 0° is considered as an empty product of zeros, its value 
is also set to 1. The following theorem holds. 
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Theorem 4 (Generalized Sylvester’s identity) Let M be a matrix of dimension 
nxm, and Mo its square leading principal submatrix of fixed order t, 0 < t < min(n, m ). 
Then, chosen s, q £ N, 1 < s < min(n — t,m — t), such that min (n,m) = t + qs, the 
following identities hold for 0 < k < q 


det Mfc 
[det M 0 f~ 1)k 


k -1 


[ [det 


i— 1 
i odd 


k -2 


n [ dets . 


—i 


l(s—1) 


1=0 
j even 


for k even, 


(15) 


fc-i 


[det Bj 


]G-i) k 


i =0 


det M k • [det Mq] 


(s—l) fc _ i even 


k—2 


for k odd. 


(16) 


II [ (let B :. 


l(s-i)' 1 


3 = 1 
j odd 


Proof. We make the proof by mathematical induction over k. 

For k = 0, (1151) is trivially satisfied since, in the right hand side, we have a ratio of 
empty products. 


For k = 1, since the product in the denominator of the right hand side is an empty 
product, we have 

det Mi • [det Mq ] s ~ 1 = detl?o, 

' 1,..., t + 

1, . . . , t + S J 


where Mi = M ( ’ ’ ’ ’ ^ + S j, detMo = af t 1 ^, and Bq is the matrix with elements 

\ 1 > ■ ■ ■, t + sj 


4? = M 


1 ,...,t,t + i 
l,...,t,t + j 


M 


= Ot+jt+j, for 1 < i,j < s. So it is exactly the Sylvester’s 

identity ([2]) applied to Mi, and (fl6l) holds. 

In the inductive step we shall prove that, if (fl6l) holds for k odd, then (fl5l) holds 
for k + 1 even, moreover if (1151) holds for k even, then (1161) holds for k + 1 odd. 


Assume that (fl6l) is satisfied for k < q odd. We apply (J2J) to the matrix Mfc + i, of 
order (m*, + s), with detM^ = aiff k k m ] j fixed, and we have 


det Mfc + i • [det Myff 1 = det-Bfc 


(17) 


where B ^ is the matrix with elements as in (1141) . Substituting in (1171) the expression 
for det Mfc obtained from (fTHl) , leads to 


det M fc+ i 
[det M 0 f~ 1)k+1 


k -1 


1 s-1 


II [det Sj] (s_1) 


i=0 
i even 


k—2 


n i detB , 


l(s-l ) k ~ 1 - j 


3 =1 
j odd 


det B k , 
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and then we get 


det M k+ 1 
[det M 0 ] (s_1) + 


which is (1151) for k + 1 even. 


3 =1 
j odd 

II [detBi]'- 1 '*-* 

i=0 
i even 


Assume, now, that (1151) holds for k < q even. By using ([2]) , the relation (ED holds 
again. Substituting in (fT71) the expression for det M k , obtained from (fl5l) . gives now 


det M k+1 • [det M 0 ] {s ~ 1)k+1 


k—1 


-i s —1 


[ [det B q 


l(s-l)' 1 


i =1 
i odd 


k -2 


[] [det Bj]( s -^ 3 


3 =0 
j even 


det B k , 


and then 


det M k+ \ ■ [det Mq] <s 


which proves (|16l) for k + 1 odd. 

□ 


n [ det ^ 




3=0 
j even 


k -1 


[det B, 




i= 1 
i odd 


So, now, we have established a relation between determinants obtained by bordering 
a leading principal submatrix of order t of a matrix M, with blocks of s rows and s 
columns, m = n (i.e. s = n — t), (fTHl) becomes exactly the classical Sylvester’s 
determinantal identity Q. 

Remark: When s = 1, we have det M k = detl?fc_i = for 0 < k < min(n — 

t,m — t), and the B k _\ are scalars. 


4.1 Applications 


Obviously, from the computational point of view, if s is large, the formulae given in the 
Theorem 0] are not very interesting. So, let M a square matrix of order n, and consider 
the particular case where s = 2, chosen t such that n — t is even, q = (n — t)/ 2, even or 
odd. We set 


M = 


^ 4 ii A12 \ 

^21 ^-22 J 


where An is a square matrix of order t, A 22 is a square matrix of order (n — t) = 2 q, 
A 12 , A 21 matrices of dimension t x (n — t) = t x 2 q. 
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Formulae (fT5l) and (fl6l) . for k = q, become 


det M 
det An 


q -1 9-2 

J det Bi / JJ det Bj , 

i=l j =0 

i odd j even 


q even, (18) 


q-l q -2 

det M - detail = det Bi / JJ det Bj , q odd, (19) 

i =0 j =1 

* even j odd 

where the Bj matrices, computed as in (1141) . are all of order 2. 

Consider the particular case where <7 = 1. Let D = An be a square matrix of order 
n — 2, 

A 22 = g 1 1 ’ ^ 12 = ( c f ) > ^21 = ( b h ) , 

with a, e, <7, l scalars, and b,c, f, h column vectors of length (n — 2). Formula (1191) 
gives 

det Mi ■ det D = det -Bo, 

and since, in this case, 


det Bq = 


D c 


D f 

b T a 


b T e 


D c 


D f 

h T g 


h T l 


we recover exactly the rule (J5|). 

Let now M' a square matrix of order n, partitioned as 

(A B t E \ 

M' = C D F I , 
\G t H t L 


( 20 ) 


where the block D is a square matrix of order t, A, E. G , L are square matrices of order 
2, and B, C, F, H are of dimension t x 2, and n = t + 4. Setting 

/ D C F \ 

M = B t A E , 

\h t g t l ) 


we have det M' = det M. So considering s = q = 2 , m gives 

det M det Bi 
det D det Bq ’ 

with Bq and Bi of order 2. 

For the above matrices we denote 

B = ( bi b 2 ),C=(ci c 2 ),E=(ei e 2 ) , F = ( fi f 2 ), 
G = ( gi g2 ), H = ( hi h 2 ) , 

A = ( ciij ), L = ( kj ), for i,j = 1,2, 















Generalizations of Sylvester’s determinantal identity 


18 


where b;, c i; e*, fj, g j, hj, i = 1,2, 

D ci 

bi T an 


are column vectors. We have 


D 

n J 


bi T 


C2 

a 12 


an 

Cl 


bi T 

D 


det Bq = 


14 ci 

b2 T 021 


D c 2 

b2 T 022 


021 b 2 T 

ci D 


012 

bi T 

C2 

D 

022 

h 2 T 

C2 

D 


detili = 


DC fi 
A ei 
hi T gf hi 

DC fi 
A ei 
h 2 T g^ Z21 


T> C f 2 
B t A e 2 
hi T gf Z 12 

D C f 2 
A e 2 
h 2 T g^ Z 22 


A B t ei 


A B t e 2 

CD fi 


CD f 2 

gi hi T Zn 


gf hi T Z i2 

A B t ei 


A e 2 

C D fi 


C D f 2 

g 2 h 2 T Z 2 i 


g-f h 2 T Z 22 


So this application is a generalization of the Case 2 of the Section [2] 


4.1.1 Examples 


Example 1: Let n = 10, t = s = 2, q = (n — t)/2 = 4 even. Thanks to (1181) . we have 

det M det B\ • det B 3 
det An det Bq ■ det B 2 

If n = 10, t = 6, s = 2, q = {n — t)/2 = 3 odd, thanks to (fl9l) . we have 

det Bq ■ det B 2 


det M ■ det An = 


det ili 


Example 2: We consider a square matrix AT of order n = 8 partitioned like in (1201) . 

r , f3,... ,6^ 


with D = M 


A = M‘ 
B t = AT 


3,...,6 

t(h 2 


Then t = 4, s = q = 2, and 

,//%2 


1,2; 5 

1,2 


E = AT 


7,8/ ’ 


, H T = M’ 


v 3,...,6 

So we can apply (12IT) . and we obtain 


det AT 


7,8 

3,...,6 


G = M 


,C = M‘ 


,(7, 

1,2; ’ 

, f 3,... ,6 

1,2 


^7 


L = Ai 
,F = M‘ 


if 3,..., 6 
7,8 


det H 


AT 

1 ,.. 

!,■■ 

- , 6,7 

- , 6,7 

• M ’ 

1 ,... ,6,8 
1 ,... ,6,8 

-AT 

'!,■■■ 

.!,■■■ 

'c- 00 1 

<0 

• AT 

1 ,..., 6,8 
1 ,..., 6 ,7 

M 1 

1 , 3 , 

...,6 

■ AT 

2 , 3 ,.. 

.,6 

-AT 

1 , 3 , 

..,6 

• AT 

2 , 3 ,.. 

• , 6 

1 , 3 , 

• • •, 6j 

2 , 3 ,. 

. ,6_J 

2 , 3 , 


1 , 3 ,. 

• ,6j 


5 Conclusions 


In this paper we presented several generalizations of the Sylvester’s determinantal iden¬ 
tity, proposed from various authors, and described here in a unified way. We started to 
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study those formulae because the classical formula is a very important tool for finding 
extrapolation algorithms mm for accelerating the convergence of scalar and vector se¬ 
quences, and we need to have at our disposal generalizations able to be used in building 
new extrapolation algorithms for vector sequences (a case more difficult) [8]. During 
our study, we had the idea for the new generalization proposed. 


Remark 

This is an updated version of a work presented in two congresses in 2008. 
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